Abstract. In this paper we present new observational ev- 
idence that supports the presence of an extra source of 
continuum emission in the binary system (3 Lyrae. New 
VLA and IRAM observations, together with published 
data from the literature and ISO archive data, allow us 
to build the Spectral Energy Distribution of the binary 
between 5 x 10 9 Hz and 5 x 10 15 Hz. The radio-millimeter 
part of the spectrum is consistent with free-free emission 
from a symbiotic-like wind associated with the primary 
. . . component and ionized by the radiation field of the hidden 
companion. Furthermore, we also consider the possibility 
that the observed radio flux originates from collimated 
. radio structures associated with the mass gaining compo- 
nent and its disk (Conical thermal jets). An extrapolation 
of this emission to the far-IR part of the spectrum indi- 
cates that in both cases the free-free contribution at these 
frequencies cannot explain the observations and that the 
observed infrared excess flux is due principally to the sec- 
ondary component and its associated disk. 
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1. Introduction 

Even if the eclipsing binary (3 Lyrae is one of the most 
studied stellar systems, its enigmatic spectroscopic and 
photometric behaviour is not fully understood. The cur- 
rent view is a non-degenerate, semi-detached interacting 
binary system in the phase of large-scale mass transfer be- 
tween components. The mass losing component is a B6- 
B8p II, while the unseen mass gaining component is prob- 
ably a BOV star, embedded in an optically thick accretion 
disk (Hubeny & Plavec 1991). The presence of a large 



plasma cloud, surrounding both components, has been in- 
ferre d from optical and UV emission lines (Batten & Sa- 
hade 1975 ; Hack et al. 1975 ) as well a s from the analysis 
of UV light curves (Kondo et al. 1994 ) . Recently the pres- 
ence of jet-like structures, probably related to the accre- 
tion disk, has been shown by the studies of Har manec et 
al. ( |l996| ) and confirmed by Hoffman et al. ( 199§|) . (3 Lyrae 
is also a well known radio source. First detected, in the 



early seventies, by Wade & Hjellming ( 1972 ) at 2.7 and 
8.1 GHz, it was then monitored, at two frequencies, by 
Gibson fll975p . The source always exhibited a thermal- 
like spectrum but Wright & Barlow ( 1975 ) pointed out 
that the observed slope of the radio spectrum (a=0.96) 
was intermediate between the slope expected for a sim- 
ple H II region and from a stellar wind, implying that the 
physics underlying the radio emission of (3 Lyrae is more 
complicated than assumed in either of these models. 

Recent high resolution MERLIN observations of 



(3 Lyrae at 4.9 GHz (Umana et al. 2000) have revealed an 
extended radio nebula around the system, whose bright- 
ness temperature (1.1 x 10 4 K) confirms the thermal origin 
of the radio emission. Such a nebula can be re-conducted 
to a massive wind associated with the cooler primary ion- 
ized by the ho tter s econdary if the two-winds model by 
Mazzali et al. ( 1992] ) is adopted. 

To assess a clear picture of the radio properties of 
(3 Lyrae, an analysis of its spectral energy distribution 
(SED) from radio to infrared appears to be necessary. In 
this paper we present multi-frequency VLA and IRAM 



observations that will be combined with data from the 
literature in order to investigate the origin of the radio 
emission of (3 Lyrae. 

In the following, assuming the most common notation, 
we will indicate as primary the mass losing component, as 
it is the most luminous in the optical region. 

2. Observations and Results 

2.1. The VLA data 

The observations were carried out using the VLA[] on May 
5, 1996, from 07:58 to 12:46 UT. We observed (3 Lyrae at 
five frequencies, namely 4.8 (C-Band), 8.4 (X-Band), 14.9 
(U-Band) and 22.0 (K-Band) and 43.0 (Q-band) GHz, 
using two independent 50 MHz bands. The observations 
were performed in compact (D) configuration, splitting 
the array in two different subarrays: the first one included 
all the antennas (13) equipped with Q-band receivers; the 
other subarray was used sequentially at each remaining 
frequency. This ensures the necessary sensitivity for detec- 
tion of weak sources. The adopted antenna configuration 
provides a typical beam size of ~ 12", ~ 7.5", ~ 4", ~ 2.5" 
and - 3" at 4.8, 8.4, 14.9, 22 and 43 GHz respectively. 

A typical observing cycle consisted of 15-20 min inte- 
gration time, preceeded and followed by a 2-min obser- 
vation of the phase calibrator. This basic sequence was 
repeated at least 3 times in order to improve the signal to 
noise ratio. 

For Q-band observations a slightly different observ- 
ing strategy was followed. Frequent checking for pointing- 
was performed since systematic errors may be a signifi- 
cant part of the primary beam at 43 GHz. Moreover, in 
order to minimise atmospheric effects on the phases, a 
much shorter observing cycle has been adopted. As phase 
calibrator, 1925+211 was chosen for the Q-band, while for 
the other frequencies we used 1850+284. The flux density 
scale was determined by observing 3C286. The 43 GHz 
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Band [GHz] Flux density [mjy] rms [mjy] U.T. 



4.8 4.27 0.06 10:40 

8.4 6.50 0.05 11:00 

14.9 9.60 0.20 10:00 

22 12.70 0.30 11:20 

43 13.90 0.50 10:00 



Band [GHz] 


Flux density [mjy] 


rms [mjy] 


U.T. 


85 


26.8 


0.7 


5:13 


115 


27.0 


4.0 


4:54 


215 


36.9 


4.5 


5:13 


240 


38.0 


3.0 


4:54 



flux of 1925+211 was determined relative to 3C286 by us- 
ing only those scans obtained at the same elevation as 
for 3C286. This ensures a careful amplitude calibration 
since at high frequencies there is a strong dependence on 
elevation of the antenna aperture efficiency and of the at- 
mospheric opacity. 

The data processing was performed using the standard 
programs of the NRAO Astronomical Image Processing 
System (AIPS). To achieve the highest possible signal to 
noise ratio, the mapping process was performed by using 
the natural weighting and the dirty map was CLEANcd 
down as close as possible to the theoretical noise. At each 
frequency, the source position and the flux density deter- 
mination were obtained by fitting a gaussian brightness 
distribution (JMFIT). To estimate the noise level in the 
maps we analyzed an area on the map, with dimensions of 
~ 15 — 20 6 syn , away from the phase center. Its consistency 
with the expected theoretical noise was always acertained. 

We detected (3 Lyrae at all the 5 frequencies as a com- 
pact, unresolved source. Our results are summarized in 
Table [j], where the radio flux density, with its associated 
rms and the time of the observation are reported for 
each different frequency. In spite of the fact that j3 Lyrae 
was one of the first stars to be detected in the early seven- 
ties, the measurements here reported are the first detec- 
tions of the system at high frequencies. 



2.2. The IRAM data 

Observations of the millimeter continuum of (3 Lyrae were 
obtained in 1996, between 15 and 16 of May, with the 
IRA Mp| Pl ateau de Bure (PdB) interferometer (Guilloteau 



et al. 1992). These observations were carried out in a com- 



pact configuration (D), providing a typical beam size of 
5" x 8" at 1.3 mm. Since each of the four operating an- 
tenna was equipped with a SIS dual- frequency receivers at 
1.3 and 3 mm, the observations were performed simultane- 
ously at 85 GHz and 215 GHz (on May 15) and at 115 GHz 
and 240 GHz (on May 16). The observations were made in 
double-sideband (DSB) mode, with each sideband of 500 
MHz separated by 3 GHz. 

A typical observation consisted of 20 minutes on source 
scan, plus several calibration measurements to assure an 



2 IRAM is supported by INSU/CNRS (France), MPG (Ger- 
many) and IGN (Spain) 



accurate calibration of the interferometer, bandpass and 
phase calibration and a high pointing accuracy. This se- 
quence was repeated several times, for a total of about 
four hours source integration time, in order to improve 
the signal to noise ratio. Phase calibration was performed 
by using 2013+370 while the flux scale was fixed by us- 
ing daily observations of the primary calibrators MWC349 
and 3C454.3. The IRAM data were reduced and mapped 
following the standard procedures using the software de- 
veloped by the Observatoire de Grenoble and IRAM 
(Grenoble Image and Line Data Analysis Software) 

Wc detected (3 Lyrae at all the 4 frequencies. The re- 
sults are summarized in Table ^, where the millimetric 
flux density, with its associated rms, is reported for each 
different frequency 



3. Spectral energy distribution 

(3 Lyrae has been observed by the ISO satellite and we re- 
trieved three scientifically validated spectra obtained with 
the Short Wavelengths Spectrograph (SWS) from Novem- 
ber 1997 to March 1998. The spectra were successively an- 
alyzed using the ISAP package. This software was specifi- 
cally written for the reduction and scientific analysis of 
ISO SWS and LWS (Long Wavelengths Spectrograph) 
Auto Analysis Results (AARs), which is the final output 
of the standard ISO pipeline processing. 

As a final step, we estimate the continuum flux in 6 
different channels by choosing areas free of spectral lines 
and evaluating the average flux by applying a box-car fil- 
ter of 2 \x width. This analysis was performed only in the 
spectrum with the best signal to noise ratio. 
(3 Lyrae is also positionally associated with the IRAS 
source 18482+3318 and, as pointed out by Friedmann et 
al. ( |1996[ ), has a spectral distribution in the IRAS bands 
typical of circumstellar material. 

Fig. [I] shows the Spectral Energy Distribution (SED) of 
(3 Lyrae built using our reported observations, referenced 
measurements from the literature as well as unpublished 
archive data. 

(3 Lyrae is reported in the literature as a periodically 
variable source, whose ephemeris are (Harmanec & Scholz 
1993| ): 

T = JD 2 408 247.966 + 12.91378S+ 3.87196 x 10" 6 ^ 2 (1) 
and questions may arise regarding whether it is possible 
to combine data obtained at different orbital phases. 
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In order to take in account the variability with the 
orbital phase in building the SED, we should consider only 
observations carried out at orbital phases close to those 
of our VLA observations, corresponding to orbital phase 
$ - 0.151 ±0.016. 

For the infrared and visible, measurements obtained 
at orbital phases between 0.14 and 0.16 and corrected for 
an interstellar extinction of -Eb-v = 0. m 065 (Abt et al., 
1962) have been used. 

For the ISO and IRAS spectral regions, information on 
possible variation with orbital phase is not available, but 
it is very probable that at longer wavelengths the contri- 
bution of the circumsystem material would be important 
and the resultant continuum will be less affected by orbital 
variation. Given the amplitude of light curves decreasing 
from the visible to the IR, proceeding in a conservative 
way, if we use the eclipse depths reported by Zeilik et al. 
( |1982| ) for the infrared as dispersion associated with each 
IRAS and ISO datum, we end up with an errorbar smaller 
than the symbols used in the plot of Fig. |]. 

Finally, we can built up the radio-mm spectrum quite 
confidentily, even if the source has been reported variable 
at these frequencies and the data were not obtained simul- 
taneously, since a VLA radio flux monitoring of (3 Lyrae 
did not reveal any variatio ns at least on a time-scale of a 
month or so (Umana et al. 2002| ). Moreover our millimet- 
ric measurements are in agr eemen t with those obtained at 
250 GHz by Altenhoff ct al. ( 1994 ) more than 15 years ago. 



Same secular changes of the radio flux density of 
(3 Lyrae are present, as reported by several authors. How- 
ever, in the following analysis of the spectral distribution 
we allow the flux density at 5 GHz to change between 2.9 
and 4.3 mJy, minimum and maximum measured 5 GHz 
flux density, retaining the spectral index determined in 
the present work. This leads to the conclusion that tak- 
ing secular changes into account does not affect our main 
results. 

4. Discussion 

The radio spectrum of (3 Lyrae can be well represented 
by a single power-law. The best fit of the cm data is 
S v oc ;,( -62±o.02)^ cons i s tent with an optically thick ther- 
mal source and very close to the canonical spectral index 
value as expected from a stellar wind. 

This result is quite different from what was claimed by 
Wright & Barlow ( 1975| ) who derived a spectral slope of 
0.96 from observations of (3 Lyrae in the early seventies 
carried out at 2.3 and 8.4 GHz. We attribute the differ- 
ence between this and our result to the large experimental 
errors associated with the old measurements. 



Jameson & King (1978) made the first attempt to 



tire system, foresees a black body spectral index and a 
turnover frequency at 84 GHz, which is not evident from 
our observations. 



4-1- Stellar wind/s 

The UV spectrum of f3 Lyrae has been known to be dom- 
inated by anomalous continuum and very strong emission 



lines (Hack et al. 1975). The observed emission features, 
which show unusually strong P-Cyg profiles, can be di- 
vided in two groups. The first group contains lines of 
highly ionized species, typical of hot-star winds, with a 
broad underlying emission. The other group includes lines 



of species at lower ionization stages (Aydin et al. 1988) 



A two- wind radiation-driven model has been proposed 



by Mazzali et al. (1992) to explain these two different 
kinds of mass loss features. The different characteristics 
suggest two distinct regions where they form, the first as- 
sociated with the wind of the B0 V gainer, the second in 
the wind of the B6-8 II loser. 

The possibility that the radio emission could be re- 
lated to the interactions between the two stellar winds is 
ruled out by the fact that the Mazzali et al. ( 1992] ) model 
outlines a system where only one of the winds is domi- 
nant. This would imply that the effects of possible wind 



interactions are negligible (Stevens 1995). Moreover, the 
observed radio spectrum is consistent with thermal emis- 
sion rather than non-thermal as observed in most colliding 
wind radio sources (Williams 1996[ ). 

Recently Umana et al. ( |200C ) resolved an extended 
radio nebula, embracing the entire binary system, by us- 
ing the MERLIN interferometer at 4.9 GHz. In this paper 
the authors suggested that it is possible to interpret the 
observed nebula, in the framework of the Mazzali et al. 



model the radio emission of j3 Lyrae by assuming a homo- 
geneous, spherically symmetric Hn region. This model, 
which assumes an emitting region embracing the en- 



(1992) hypothesis, as a result of the the stellar wind as- 
sociated with the BII component ionized by the hotter 
companion. Such a phenomenon is observed in Symbiotic 
systems where the associated radio nebula is the material 
ejected by one of the components which is ionized by the 
secondary star. 

The radio properties of symbiotics are well summarized 
in a scries of papers (see for example, Seaquist & Taylor 
1990). In particular, the observed radio spectra have been 
modelled in terms of a binary model, where the morphol- 
ogy of the radio emitting region is a function of the physi- 
cal characteristics of the system such as binary separation 
(a), Lyman continuum luminosity of the hot component 
(I/ uv ), mass loss rate (M) and velocity (v) of the dominant 
wind (Taylor & Seaquist |l984| TS). 

If the wind is completely ionized, the ionized region 
shows the same radio properties as a spherically symmet- 
ric ionized wind. Thus, in the framework of the TS model, 
a spectral index close to the canonical 0.6, as we observed 
in (3 Lyrae, would imply that at least a large fraction of 
the wind is ionized. A large ionization region can be writ- 
ten in the TS terminology as q Q a << a, where q Q is the 
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Fig.l 

datum from Altenhoff et al 



The SED of (3 Lyrae. Dots, diamonds and the crosses represent VLA and IRAM data from this paper plus a 
(1994). Triangles and squares indicate IRAS and ISO data, while asterisks are VJKLM 



magnitudes from Viotti et al. (1978), V magnitude has been corrected as in the text. Data points have been corrected for 
interstellar extinction following Mathis ( |1990[ ) and assuming a £"b-v = 0. m 065. The dotted line is the contribution at 
high frequencies of the free-free plus bound-free emission from the stellar wind. The stellar contribution is represented 
by an ATLAS9 (Kurucz 1993) model, with fluxes normalized to the flux in the V band (solid line). 



distance, measured in terms of the binary separation a, 
from the center of the mass-losing star and the boundary 
of the ionized region. Using this assumption, from equa- 
tion 8 of TS, we can relate the turnover frequency to the 
physical parameters of the system by the equation: 



t^-1.35 ji> 

vl 1 » 8 x 10 66 (^^)(— ) 2 

1 y a 3 /v v ' 



(2) 



where a is in cm, M in MQyr , v in km s . 
Assuming a = 58.5i? Q (Bisikalo et al. |2000D , the M 

and v values of the dominant wind (BII component, see 
Table 3) and T c va lue obtained by MERLIN observations 
(Umana etal. |2000D , we derive a turnover frequency higher 
than 760 GHz. 



The observed spectrum, which is optically thick up to 
~ 300 GHz, is thus consistent with the symbiotic inter- 
pretation. 

The far IR part of the spectrum is more difficult to 
interpret mostly because at these frequencies the contri- 
bution of the photosperic emission becomes important. 

The primary star in (3 Lyrae is spectroscopically clearly 
visible and has been classified as a B8-6 II (Balachandran 



et al. 1986). On the contrary very little is still known about 



the nature of the secondary nother than the fact that it 
should be embedded in an accreti on dis k. However, despite 
of th e numerous eff orts ( Wilson [1974 ; Linnell & Hubeny 
1996| ; Linnell et al. |l998t Linnell [2000] ), there is no stan- 
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excess flux 




Table 3. Binary and Nebula Characteristics 



Binary 



Gainer 



Loser 



T cB [K] 


32 000 


13 300 




4.68 x 10 -8 


7.16 x 10 -7 


Voofkm s _1 ] 


1470 


388 


Radio Nebula 4, 






T e [K] 


11 000 ±700 




Size (FWHM) [mas] 


(145 ±12) x (100 ±8) 





Frequency (Hz) 



Fig. 2. The excess flux of Lyrae. The line is the fit of 3) Mazzali et al. (|1992|) 
the contribution at high frequencies of star±disk system. 4) Umana et al. fl2000| ) 



dard model of an accretion disk able to reproduce syn- 
thetic light curves that fit the observed ones from IR to 
UV. 

Therefore, we can compute the expected far-IR flux 
only for the primary component by assuming, between 
1.5 x 10 16 and 1.5 x 10 12 Hz a Kurucz ( |1993| ) model with 
T e g = 13 000 K, log(<?) = 3, as appropriate for the B8-6 II 
component. 

Since the observed continuum is the sum of the pri- 
mary and secondary continuum, to normalize the atmo- 
spheric model we correc t the observed V magnitude, as 
observed by Viotti et al. (1978), by a factor 1.10. This has 
been calculated by assuming a flux ratio between the two 
components at maximum light of 3 (Harmancc [1990 ) and 
a further decrement of about 12% of the primary's flux 
due to the fact that the observations took place at orbital 
phase $ = 0.15. 

The radio emission from a stellar wind comes from 
regions quite far from the star, where the wind velocity 
has reached its terminal velocity (voo). Because of the de- 
pendence of thermal opacity on wavelengths (n\ oc A 2 ), 
emission at higher frequencies would come from regions 
much closer to the stellar photosphere, where some accel- 
eration takes place. Thus to evaluate the contribution of 
such a stellar wind to the far-IR part of the spectrum, the 
effects of the velocity law of the velocity of the mass-flow 
must be taken into account. 



Lamers & Waters (1984) have computed the energy 
distribution for a grid of expanding, isothermal stellar 
wind models. In their work they considered all the possible 
ways the extended envelope can affect the energy distri- 
bution of the central star, modelling quite accurately the 
far-IR part of the spectrum. 



In adapting the Lamers & Waters (1984) model to the 



two winds scenario of (5 Lyrae we will assume that the 
secondary component plays only the role of providing the 
radiation field necessary to ionize the wind associated with 



the primary. Following them, we thus derive the emission 
due to the star±wind system, by assuming as physical 
parameters of the stellar wind those of the dominant wind, 
tha t of t he B8-6 II component, as derived by Umana et 
al. ( E000| ). For the wind, we assume that the velocity is a 
function of the radial distance, in units of 



^ =0.01 + 0.99(1 --) 

Voo & 



(3) 



as sugg ested by spectroscopic UV observations (Castor & 
Lamers 1979). 

The obtained energy distribution, normalized at 
5 GHz, is shown in Fig. [I] (dotted line), where it is ev- 
ident that the observed flux starts to deviate from the 
expected oc ;/ 2 photospheric behavior at ~ 5 10 Hz and 
there is an extra far-IR excess flux that the presence of a 
stellar wind cannot explain. 

This excess, extracted and plotted in Fig. g, can be fitted 
by a power law: 

S v oc V l - 10±0M (4) 
in agreement, surprising, with the results from Viotti et 
al. ( 1978| ), which, however, were based only on JKLM ob- 
servations reported in the plot. As already pointed out 
by Viotti and co-authors, this result is consistent with a 
secondary hidden by an accretion disk, whose dimension 
increases with wavelength. 



4 ■ 2. Thermal jets 

The two-wind model of Mazzali has been recently criti- 
cized by Harmanec (2002), who pointed out that there is 
no evidence of wind associated with the secondary compo- 
nent, at least from emission or P Cygni lines as observed 
in the ultraviolet. Moreover, the elongated radio nebula, 
as observed by MERLIN, appears to be aligned with the 
reported optical jets, suggesting the possibility that the ra- 
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dio emission is due to thermal radio jets emanating from 
the accretion region around the hidden secondary. 

Several authors have investigated the effects on radio 
flux, spectral index and mass-loss derivation due to de- 
viations from sphericity in a extended stellar envelopes 
(Schmid-Burkg [l982| ; Reynolds |l986| ). While collimated 
thermal jets can produce radio spectra significantly differ- 
ent from that of a canonical wind, there are no spectral 
clues that can help in discriminating between a spherically 
symmetric wind or biconical constant velocity flows (wider 
thermal jets). Therefore, we have a piece of observational 
evidence that if the origin of radio emission is a thermal 
jet it must be not collimated. 

We can derive the jet opening angle (#o) from the ratio 
of the minor (9 m i n ) and major axis (9 ma x) of the radio 
nebula as derived from MERLIN measurements (table 3). 



do = 2tg~ 



69° 



(5) 



To determine the mass-loss of a thermal jet with open- 
ing angle 9q, the formula, derived in the hypothesis of 
symmetry in the wind, has to be corrected by a factor of 
?y = 0.2 9o(sini)~ 1 / 4 (eq. 21 of Reynolds, 1986), where i is 
the angle of inclination of the jet as seen by the observer 
(i ~ 85° — 88°) and 9q is expressed in radians. 
Therefore, we may derive the mass-loss from the relation: 
M' sphe = 6.7 x m-^S^iD^liv x g s )- Q - 5 M©^" 1 ^) 
where full ionization and cosmic abundances have been 
assumed, Voo is the terminal velocity of the wind, Se cm the 
observed flux density at 6 cm in mJy, Z?kpc the distance 
of the system, in kpc. gg represents the free-free Gaunt 
factor that, following Leitherer & Robert (1991), can be 
approximated with: 

T 3/2 

gs = 9.77(1 + 0.13 log ) (7) 

v 

where T, in Kelvin, is the wind temperature. 
By assuming a typical stellar wind velocity ofa hot B star 
for the secondary component (v ~ 1500 km s _1 ) and con- 
stant temperature equal to the brightness temperature, 
derived by radio observations (T — Tb), we obtain a mass 
loss rate of 

M' sphe = 2.0 x 10- 6 Moyr- 1 . (8) 
which, corrected for thermal jets yields a mass-loss of 



jet 



7]M sphe = 5.2 x 1Q- ' M yr" 



(9) 



A wide thermal jet will emit a radio flux whose frequency 
dependence is indistinguishable from that emerging from 
a spherically symmetric stellar wind. Therfore, we can ap- 
ply the model of Lamers & Waters (1984), introducted 
in the previous section, to derive the emission due to the 
star+jets system where, this time, the physical parameter 
of the stellar wind associated with the gainer (see table 
3) and mass-loss corrected for thermal jets have been as- 
sumed. 

The obtained energy distribution, normalized at 5 GHz, 
is showed in Fig. || (dotted- line) , where the contribution 
of the primary (as indicated by the KURUCZ model) has 
been also taken into account. For a better comparison, 



the model, obtained in the case of symbiotic winds (see 
section 4.1) is shown in the figure as a dashed line. The 
model stops at v ~ 10 14 Hz as this is the highest frequency 
of the Lamers & Waters (1984) tabulation. 

As already noted for the symbiotic wind hypothesis, 
there is an extra far-IR excess flux that a thermal jet, 
associated with the secondary, is not able to explain. The 
extracted excess flux is well fitted by a power-law, with 
a spectral index a = 0.91 ± 0.06, still consistent with the 
contribution of an accretion disk, whose density decreases 
with radius. 



5. Conclusions 

j3 Lyrae is one of the best studied stellar systems: still no 
reliable model of the system, able to reproduce the huge 
amount of observational data, is available. 
Very recently, a new generation of models, that take into 
consideration the details of the physics of the accretions 
dis k, hav e been presented (Linncll ct al. 1998 : Bisikalo ct 
al. 2000 ). Even if they provide a good model for the light 
curves in the near and far-UV and in the visible, they 
still do not predict the shape of the infrared light curve, 
and, in particular, they do not reproduce the secondary 
minimum that becomes deeper than the primary mini- 
mum at A > 4.8/xm. This leads to the conclusion that an- 
other component, besides stars plus accretion disk, must 
contribute to the observed flux, and Linnell (2000) pro- 



poses that the extra source of continuum radiation may 
be Thomson scattering of radiation from the gainer. 

In this paper we have presented new observational evi- 
dence that support the presence of an extra-component in 
the already quite complex binary system. This extra com- 
ponent, which is the origin of the observed thermal radio 
flux, can be due to an extended stellar wind ionized by 
the strong ultraviolet flux of the secondary component or 
related to collimated structures associated with the gainer 
(conical thermal jets). 

We further evaluated the contribution of this compo- 
nent in the far-IR to extract the spectral energy distribu- 
tion of the stars plus accretion disk system. Our results, 
in both hypotheses, show a power-law distribution up to 
IRAS frequencies indicating an accretion disk with a non- 
uniform density distribution and whose size varies with 
wavelength. To quantify this, the collection of new, good 
quality infrared light curves would be highly valuable. 

The stellar wind or thermal jets, as appearing from 
the obtained spectrum and its modelling, is optically thick 
up to IRAS-ISO frequencies. Still, at those frequencies we 
measure a significant excess that we attribute to the disk 
since it is in agreement with the near-IR data. We may ask 
what kind of effect the radio nebula may have on the ob- 
servability of the inner disk. In both the cases, the size of 
the radio source depends on frequency as i? ra dio oc j/~ - 7 . 
Thus, scaling the size measured at 4.8 GHz with MERLIN 
(see Table 3) we obtain that the radio nebula reaches di- 




Fig. 3. Same as Fig. [T]. The dotted line is the contribution at high frequencies of the free-free plus b ound-free emission 
from a biconical outflow plus the secondary contribution represented by an ATLAS9 (Kurucz 1993| ) model (continuos 
line), with fluxes normalized to the flux in the V band (solid line). For comparison the fit in the case of symbiotic 
wind is also shown (dashed line) 



mensions comparable to those of the entire binary system 
at frequencies of the order of 10 13 Hz, close to the IRAS 
spectral region. 

This provides important constraints to the radio source 
associated with the binary system, whose morphology 
should allow one to observe also the disk. The other pos- 
sibility, that the disk is surrounding the nebula, can be 
discarded, since a clear eclipse is observed in the infrared. 

Finally, we would like to point out that even in the case 
of stellar wind, the fact that the radiation field necessary 
to ionize the wind is provided by the secondary compo- 
nent, which is embedded in the thick accretion disk (Lin- 
ncll 2000), suggests that the ionization would take place 



mainly in the polar regions. Thus, in both the cases, it is 
quite probable that the nebula has a bipolar morphology 



that the MERLIN observations (Umana et al. 2000) were 



not able to resolve. This bipolar morphology, which should 
have a wide opening angle as suggested by the observed 



spectral index (Rodriguez et al. 1990), is compatible with 
the possibility of observing the inner accretion disk. 
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